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^ ■ Abstract 

We study the antiferrodistortive instability and its interaction with ferroelec- 
^^ . tricity in cubic perovskite compounds. Our first-principles calculations show 

^^ . that coexistence of both instabilities is very common. We develop a first- 

o 

principles scheme to study the thermodynamics of these compounds when 
Q^ , both instabilities are present, and apply it to SrTiOs. We find that increased 



:a 



X 



pressure enhances the antiferrodistortive instability while suppressing the fer- 
roelectric one. Moreover, the presence of one instability tends to suppress the 
other. A very rich P-T phase diagram results. 
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The great fascination of the cubic perovskite structure is that it can readily display a 
variety of structural phase transitions, ranging from non-polar antiferrodistortive (AFD) 
to ferroelectric (FE) and antiferroelectric in nature |1|]. The competition between these 
different instabilities evidently plays itself out in a variety of ways, depending on the chemical 
species involved, leading to the unusual variety and richness of the observed structural phase 
diagrams. Moreover, all the phase transitions involve only small distortions from the ideal 
cubic structure, and are therefore appealing objects for experimental and theoretical study. 
However, our microscopic understanding of the chemical origins of these instabilities and of 
their interactions is still very limited. 

Thus, there is a pressing need for accurate, chemically specific investigations of the 
structural energetics of these compounds, leading to a detailed understanding of the phase 
transition behavior. Previous phenomenological model Hamiltonian approaches pj-Q have 
been limited by oversimplification and ambiguities in interpretation of experiment, while 
empirical and nonempirical pair-potential methods [0 have not offered high enough ac- 
curacy. First-principles density-functional calculations have been shown to provide accurate 
total-energy surfaces for perovskites as regards FE distortions P^TI]|. However, to our knowl- 



edge, there have been no previous first-principles studies of AFD distortions, and therefore 
no detailed microscopic theories of the phase transformation behavior. 

Here, we build upon previous work in which a fully first-principles scheme was used to 
study the FE transitions in BaTiOa, leading to an accurate microscopic understanding of 
the phase transition sequence [|lT|. In the present work, we develop a similar approach which 
is capable of treating simultaneously the FE and AFD degrees of freedom, allowing for the 
first time a detailed ab-initio study of the phase behavior for perovskites in which both 
instabilities are present. We present systematic calculations of the susceptibility against 
R-point zone-boundary AFD modes for a set of eight compounds, demonstrating that the 
AFD instability is very common. Then, we briefly describe our first-principles scheme for 
studying finite-temperature properties, and apply it to SrTiOs. We study the evolution of 
the phonon instabilities with temperature, and calculate the P-T phase diagram. In so 



doing, we compute the interactions between the AFD and FE instabihties, and expose their 
imphcations for the thermodynamic properties. 

The high-symmetry ABO3 perovskite structure is simple cubic with O atoms at the face 
centers and metal atoms A and B at the cube corner and body center, respectively. The 
two most common instabilities result from a softening of either a zone-center polar phonon 
mode (FE), or a non-polar zone-boundary mode (AFD) involving rigid rotations of oxygen 
octahedra. These modes are illustrated in the left and right insets, respectively, in Fig. |l|. 
BaTiOa is a classical example of the first type, while the best-known example of the second 
kind is the T = 1057^ transition in SrTiOa |^, which results from a softening of a r25 phonon 
at R [(lll)7r/a]. 

The stability of perovskite compounds against -R-point phonon distortions can be ex- 
pressed in terms of a stiffness k^ = {l/2)d'^E/d(j)'^, where is the rotation angle of the 
oxygen octahedra. To obtain k,^, we perform frozen phonon calculations using density- 
functional theory within the local-density approximation (LDA) and Vanderbilt ultra-soft 
pseudopotentials [jI2[. In Table |, we list values of k^ for a set of eight compounds, calcu- 



lated at the experimental lattice constants |T^ as listed in Ref. |]TD|]. Negative values indicate 
instability to -R-point phonon distortions. 

Table | shows that the tendency towards AFD instability is strongly correlated with 
trends in ionic radii. Such trends in an ABO3 compound are conventionally described by a 
tolerance factor t = (r^ + ro)/V^{rB + ^o)- Values for t are given in Table | using the ionic 



radii of Ref. Jl^. We find that kr is almost monotonic with t, i.e., a larger A or a smaller 
B atom tends to stabilize the cubic structure. This simple behavior contrasts with the case 
of the ferroelectric instability, where covalent interactions play an important role |]15| . 



Inspecting Table |, we see that the two compounds BaTiOa and KNbOa are clearly stable 
with respect to AFD distortions, consistent with experimental observations. (Both materials 
undergo a similar series of FE transitions.) On the other extreme, we find that CaTiOs, 
PbZrOs, and NaNbOs have a strong AFD instability. All three compounds are also predicted 



to have FE instabilities [|10l, consistent with the observation of complex phase diagrams 



and high transition temperatures in all three cases. Finally, our calculations for SrTiOs, 
PbTiOa, and BaZrOs show a weak AFD instability. PbTiOs is observed to go through 
a weak unidentified transition at T=180K |T^ which could be AFD related. BaZrOs is 
observed to remain cubic down to T=0; the weak instability predicted by our calculation 
could be suppressed by quantum zero-point fiuctuations. For SrTiOs, we predict a weak 
AFD instability consistent with a low Tc of 105K observed for its cubic-to-AFD transition. 

The above calculations indicate that coexistence of FE and AFD instabilities is very 
common in perovskites. To study the consequence of such a situation, we have chosen to 
study the case of SrTiOs in depth. Our first-principles scheme can be explained briefiy as 
follows. The energy is Taylor-expanded in low-energy distortions, with expansion parameters 
determined from LDA calculations. The resulting Hamiltonian is studied using Monte Carlo 
(MC) simulations. The low-energy distortions we included are those connected with zone- 
center F15-like modes, zone-boundary AFD-like modes, and strain. To do this we construct 
a FE "local mode" such that a uniform arrangement of local mode amplitudes f; reproduces 
the softest zone-center Fis (FE) modes (/ is a cell index). Similarly, we construct an AFD 
local mode (a local rotation of an oxygen octahedron) so that a staggered arrangement of 
amplitudes a; reproduces the r25(-R) mode. Finally, the local strains are represented in 
terms of a displacement vector u;. 

Thus, we have three vector degrees of freedom f^, a;, and U; per cell. The energy terms 
retained in our Taylor expansion of the potential energy are: (i) on-site self energy, up to 
quartic anharmonic order for f; and a;, and up to harmonic order for u; (elastic energy); (ii) 
harmonic intersite interactions between f; (including long-range dipole-dipole interactions) 
and between a; (short-range only); and (iii) on-site coupling energy to the lowest order 
between a^ and u;, between f; and u;, and between f; and a^. The determination of the 
expansion parameters involves LDA calculations for supercells containing up to 20 atoms 
with low symmetry, using ultra-soft pseudopotentials |12[ . The details of the Hamiltonian, 
the first-principles calculations, and the values of the expansion parameters will be presented 
elsewhere fl^ . 



To obtain the structural and thermodynamic properties, we perform MC simulations 



on an L X L X L cubic lattice with periodic boundary conditions |]T8|. The identification 
of different phases can be made by monitoring the FE order parameter f (F) (the Fourier 
transform of f^ at k = F), and similarly the AFD order parameters a(/2) and a(M) [M = 
(110)7r/a]. a(M) is found to remain small for SrTiOs, and will not be discussed further. 

We first investigate the ground-state structure for SrTiOs as a function of hydrostatic 
pressure. We find it convenient to run the MC calculations at L = 4 at T=0.1K (finite- 
size and hysteresis effects are not important at low T). The calculated order parameters 
a(i?) and f (F) are shown in Fig. |l]. Zero pressure in the figure corresponds to the LDA- 
calculated equilibrium lattice constant, which is about 1% too small. Since both the FE and 
AFD instabilities are sensitive to lattice constant, comparison with the experimental phase 
diagram is best made with the zero of the pressure axis shifted by Pq=—5A Gpa (see arrow in 
Fig. 0), the value which restores the experimental lattice constant. From Fig. ||, we see that 
pressure has opposite effects on a(_R) and f (F), and that as a function of pressure, the ground 
state of SrTiOa can have four phases. The cubic phase, which is stable at high temperature, 
is not present. At high pressure, only one component of a(i?) is non-zero, indicating an AFD 
tetragonal structure {14/ mem). As P is lowered, the corresponding (z) component of f(r) 
becomes non-zero, and the structure transforms to tetragonal with FE and AFD (Mem,). A 
further decrease of pressure creates a low-symmetry monoclinic structure (Pb), in which all 
components of f (F) and a(i?) are non-zero. Finally, below —8 Gpa the structure becomes 
FE rhombohedral {Rim). We see that the coexistence of zone-center and zone-boundary 
instabilities creates many different phases and complicated structures, even at T=0. 

At finite temperature, the behavior becomes even more interesting. We first show our MC 
simulation for P = Pq (— 5.4Gpa) and L = 12. We start at high temperature and decrease 
T in small steps, allowing the system to reach equilibrium at each step. The hysteresis and 
finite-size effects on the transition temperatures are found to be negligible. In Fig. 0, we 
show the order parameters f (F) and Si{R) as a function of T. (Since the order parameter 
vectors may rotate, what we actually show are the averaged maximum, intermediate, and 



minimuni components of each vector.) Naturally, the system is found to adopt the cubic 
structure at high temperature. As T is reduced, a transition to an AFD tetragonal structure 
occurs at 130K, as indicated by a strong increase of az{R). A second transition occurs at 
T=70K to a FE tetragonal structure, below which /^(r) > 0. At very low temperature 
(lOK), the system transforms to the low-symmetry monoclinic structure. 

Comparing with experiment, we see that our cubic-to-FE(T) transition at 130K corre- 
sponds very well to the observed one at 105K |]l9l . Our observations of additional transitions 



to AFD-I-FE phases at 70K and lOK are not, however, in direct accord with experiment. 
Instead, they agree with the observed softening of the FE polar phonons, which would ex- 
trapolate to a FE transition close to 40K pO|] or 20K |^. It has been speculated that the 



absence of a true FE phase at T=0 is a result of quantum fluctuations of atomic positions, 
leading to crossover into a "quantum paraelectric phase" at very low temperature p0|-p2 



Our inability to obtain agreement between the classical MC theory and experiment at T=0 
lends additional support to this conclusion. 

To construct a P-T phase diagram, we have carried out a series of similar cooling-down 
simulations at different pressures. As shown in Fig. |^, there are at least seven different 
phases present. At strong negative pressure, SrTiOs behaves rather like BaTiOs, with a 
cubic — i> tetragonal -^ orthorhombic -^ rhombohedral sequence of transitions on cooling. 
Increasing the hydrostatic pressure tends to stabilize the AFD state and destabilize the FE 
one. The pressure coefficient dT^/dP = 28K/Gpa at Pq agrees well with the experimental 
value of 25K/Gpa p3[. At very high pressure, the system undergoes a single transition 
to a tetragonal AFD structure. In the intermediate regime, the presence of both kinds of 
instability creates a variety of phases, including the complicated monoclinic structure. The 
ordering of the FE and AFD transition temperatures reverses ~1.5 Gpa below Pq (hashed 
area in Fig. ^. In this critical region the AFD and FE transition temperature change 
dramatically, and the system may possess some interesting characteristics (e.g., extreme 
dielectric properties). 

The dramatic reversal of the AFD and FE transition temperatures in the hashed region 
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of the P-T phase diagram suggests the presence of a competition between the two insta- 
bihties. Our first-principles theory confirms this and provides microscopic insight into the 
competition. The FE and AFD instabihties affect each other mainly through the on-site 
anharmonic coupling, and through their mutual coupling to the elasticity. In SrTiOs, the 
on-site coupling is found to lead the FE and AFD modes to suppress one another, while the 
coupling through strain tends to stabilize tetragonal phases relative to other phases. Our 
calculations show that the former effect dominates. 

One way of quantifying the importance of this competition is to compare with what 
would happen if the FE or AFD degrees of freedom were artificially frozen out. We find 
that at Pq, the AFD transition temperature would be 25% higher if all f; were frozen to 
zero; conversely, the FE C-T transition would be 20% higher if all aii=0. At T=0, freezing 
f/=0 reduces the cubic-to-AFD transition pressure from —8 to —11.8 Gpa, while freezing 
a/=0 increases the cubic-to-FE transition from —1.5 to 0.8 Gpa. Thus, we see clearly that 
the FE and AFD instabilities compete with and tend to suppress one another. Because of 
this competition the T(A,F) phase at Pq is only slightly more stable than the T(A) phase, 
even at T=0; the energy surface relative to the FE distortion takes the form of a very long 
and shallow double well. This may help explain the observed suppression of the ferroelectric 
phase by quantum fluctuations ||20|^2|. 

Much of the interesting portion of our phase diagram appears to the left of Pq, i.e, at 
negative (inaccessible) physical pressures. It would be interesting, therefore, to consider 
compounds such as CaTiOs or NaNbOs which are FE at Pq, and study AFD instabilities at 
elevated P. While the exact details of our phase diagram for SrTiOa should not be expected 
to carry over to other perovskites, we expect the general features to persist, especially 
the tendency of the FE and AFD instabilities to suppress each other and the presence of 
complicated phase diagrams with numerous phases. 

In conclusion, we have performed a fully first-principles study of the finite-temperature 
properties of perovskite compounds with both FE and AFD type instabilities. We find that 
AFD instabilities are almost as common as FE ones in cubic perovskite compounds. For 



SrTiOs, our calculated P-T phase diagram shows that the FE and AFD instabilities have 
opposite trends with pressure. The anharmonic on-site coupling between order parameters 
causes the AFD and FE instabilities to tend to suppress one another. 

We thank R.D. King-Smith and K.M. Rabe for discussions. This work was supported 
by ONR grant N00014-91-J-1184. Partial supercomputing support was provided by NCSA 
grant DMR920003N. 
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TABLES 
TABLE L Calculated stiffness k of R-point AFD phonon mode (in Hartree), and tolerance 
factor t, for several perovskite compounds. 





K^ 


t 




K^ 


t 


BaTiOa 


0.295 


1.07 


SrTiOs 


-0.042 


1.01 


KNbOs 


0.242 


1.06 


NaNbOg 


-0.133 


0.97 


BaZrOa 


-0.021 


1.01 


PbZrOa 


-0.324 


0.97 


PbTiOg 


-0.037 


1.03 


CaTiOa 


-0.375 


0.97 
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FIGURES 
FIG. 1. T=0 order parameters vs. pressure for SrTiOs. Solid and dashed lines denote Cartesian 

components of f(r) and a{R) respectively. Phases are labeled by lattice symmetry (R=rhom- 

bohedral, M=monoclinic, T=tetragonal, 0=orthorhombic) and by instabilities present (A=anti- 

ferrodistortive, F=ferroelectric). Dotted lines denote phase boundaries. Vertical arrow indicates 

theoretical pressure Pq at which the lattice constant matches the experimental P=0 one. Left 

inset: sketch of displacements leading to R(F) phase (Sr is omitted for clarity). Right inset: same 

for T(A) phase. 

FIG. 2. Order parameters of SrTiOa vs. T at Pq = —5.4 Gpa. Upper panel: averaged largest, 
middle, and smallest Cartesian components of a{R). Lower panel: corresponding quantities for 
f(r). Phase labels are the same as in Fig. 1. 

FIG. 3. P-T phase diagram of SrTiOs. Hashing indicates the critical region where dramatic 
changes occur. Vertical dash-dotted line indicates the pressure Pq corresponding to experimental 
P=0. Phase labels are the same as in Fig. 1. 
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